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Abstract 

The electromagnetic interactions of a relativistic two-body bound state are 
formulated in three dimensions using an equal-time (ET) formalism. This 
involves a systematic reduction of four-dimensional dynamics to a three- 
dimensional form by integrating out the time components of relative momenta. 
A conserved electromagnetic current is developed for the ET formalism. It is 
shown that consistent truncations of the electromagnetic current and the AW 
interaction kernel may be made, order-by-order in the coupling constants, such 
that appropriate Ward-Takahashi identities are satisfied. A meson-exchange 
model of the NN interaction is used to calculate deuteron vertex functions. 
Calculations of electromagnetic form factors for elastic scattering of elec- 
trons by deuterium are performed using an impulse- approximation current. 
Negative-energy components of the deuteron's vertex function and retarda- 
tion effects in the meson-exchange interaction are found to have only minor 
effects on the deuteron form factors. 



I. INTRODUCTION AND SUMMARY OF RESULTS 

Experiments being performed at the Thomas Jefferson National Accelerator Facility 
(TJNAF) are designed to test our understanding of deuteron properties at space-like mo- 
mentum transfers comparable to the nucleon mass. In building theoretical models of these 
processes, relativistic kinematics and dynamics would seem to be called for. Hence, con- 
siderable theoretical effort has been invested in constructing relativistic formalisms for the 
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two-nucleon bound state that are based on an effective quantum-field-theory Lagrangian. 
If the usual hadronic degrees of freedom appear in the Lagrangian, then this strategy can 
be used to obtain a logical extension of the standard non-relativistic treatment of the two- 
nucleon system. The central goal of such an approach is to capture the relativistic effects 
in the two-nucleon system. At the same time, for applications to electromagnetic physics, 
it is critical to embed the Ward-Takahashi identities that imply current conservation for 
electromagnetic interactions in the theory. 

This paper focuses on a three-dimensional (3D) formalism that, in principle, is equivalent 
to the four- dimensional (4D) Bethe-Salpeter formalism. This approach has been developed 
in two recent papers [0,0. We review the formalism for relativistic bound states, and pro- 
vide an extension that ensures nonsingular behavior of the interaction in frames where the 
total momentum of the bound system is nonzero. We also construct the corresponding 
electromagnetic current. Calculations of elastic electron-deuteron scattering are performed 
based upon the impulse approximation and the results for the observables A, B and T 20 are 
presented. 

The Bethe-Salpeter equation, 

T = K + KG T, (1.1) 
for the four-dimensional NN amplitude T provides a theoretical description of the deuteron 



that includes relativity. In Eq. (IT), K is the Bethe-Salpeter kernel, and Go is the free two- 
nucleon propagator. In a strict quantum-field-theory treatment, the kernel K includes the 
infinite set of two-particle irreducible NN — ► NN Feynman graphs. Numerical methods to 
calculate the 4D t-matrix of Eq. fll.lj ), with an infinite kernel containing all crossed ladders, 
have been developed and demonstrated for a simple scalar field theory by Nieuwenhuis and 
Tjonl. 

For the two-nucleon system, such an application of the full effective quantum field theory 
of nucleons and mesons is not only impractical, it also is inappropriate because hadronic 
degrees of freedom are not fundamental. Rather, the Bethe-Salpeter formalism serves as 
a theoretical framework within which a relativistic effective interaction may be developed. 
This is entirely analogous to the way that the Schrodinger equation serves as a framework 
for development of a non-relativistic potential that describes the NN phase shifts. Since 
the NN interaction is an effective one, it is equally appropriate to develop the relativistic 
effective interaction within an equivalent three-dimensional formalism that is obtained from 
the four-dimensional Bethe-Salpeter formalism via a systematic reduction technique. 

In Ref. [[IJ, a 3D reduction of the Bethe-Salpeter formalism was developed such that the 
resulting equations involve the same propagator as appears in the Salpeter equation [|J, 

T 1 = Ki + K 1 (G Q )T x . (1.2) 

The three-dimensional propagator (Go) is obtained by integrating over the time-component 
of relative momentum, 

(Go)= J^G (p;P). (1.3) 
The 3D kernel is defined by solving two coupled equations, 
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K X = {G Q )- X {G Q KQ)(G )- X , 



(1.4) 



which is three-dimensional, and 

Q = G + G {K-K 1 )g, (1.5) 



which is four dimensional. In fact, the condition ( |1.4|) is obtained by demanding that 
(g) = (Gq). (See Ref. |1[] for details. Other works which have considered a similar formalism 
include Refs. |0-§|.) This formalism is systematic in the sense that, given a perturbative 
expansion for the 4D kernel, K, a perturbative expansion for the 3D kernel, K±, can be 
developed. However, it is necessary for (Go) to be invertible in order for the 3D reduction 
to contain the full content of the 4D theory. 

As mentioned above, when Go is the standard two-particle propagator of the Bethe- 
Salpeter formalism, (Go) is the 3D Salpeter propagator. For spin-1/2 particles it is, 

A+A 2 + AfA; 
(G ° > = P°- ei -e 2 -p» + fl+e2 ' (L6) 

where A ± are related to projection operators onto positive and negative-energy states of the 
Dirac equation, P° is the total energy of the two-body system, and 

e l = (tf + m*) 1/2 . (1-7) 

Here, and throughout the paper, these single-particle energies are to be understood as having 
an infinitesimal negative imaginary part. This defines our ie prescriptions. Note that for 
spin-half particles, the equal-time propagator of the Salpeter equation is defined only on 
two of the four sectors of the Dirac space of two particles. It is not invertible. Although 

consistent equations can be projected out for the ++ and components, one must set 

the H — and — h components to zero for consistency. Consequently, the 3D reduction does 
not have the full content of the 4D theory. If the 4D and 3D theories are to have the same 
dynamical content then we must include graphs involving H — states, such as the the time- 
ordered Z-graph of Fig. [I], in both. This graph is not contained in the ladder Bethe-Salpeter 
equation scattering series. 

The absence of this mechanism from the ladder Bethe-Salpeter equation is related to the 
fact that if we allow one particle's mass to tend to infinity, the propagator does not reduce 
to the Dirac propagator for the other particle. Consequently, Eq. (|1 . 2|) does not possess the 
correct one-body limit if any finite set of graphs is chosen for K\. 

If a scattering equation with a kernel which contains only a finite number of graphs is 
to possess the correct one-body limit, two distinct criteria must be satisfied. First the 3D 
propagator should limit to the one-body propagator for one particle (either Dirac or Klein- 
Gordon, depending upon the spin) as the other particle's mass tends to infinity. Second, 
as either particle's mass tends to infinity, the equation should become equivalent to one in 
which the interaction, K±, is static. 

In fact, Eq. (|1.2| )'s lack of either of these properties springs directly from the Bethe- 
Salpeter equation (|1 . 1|) not having the correct one-body limit if the kernel does not include 
the infinite set of crossed-ladder graphs. As alluded to above, solution of the Bethe-Salpeter 
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equation with such a kernel is impractical in the NN system. Nevertheless, the desired prop- 
erties may be obtained by using a 4D integral equation for K to reorganize the contributions 
to the kernel of the Bethe-Salpeter equation as follows, 

K = U + UG C K. (1.8) 

Given a suitable choice for Gc, this equation defines a reduced kernel, U, in terms of the 
original kernel, K. The propagator, Gc, is chosen so as to separate the parts of the kernel, 
K, that are necessary to obtain the one-body limits from the parts that are not. Once this 
is done, U may be truncated at any desired order without losing the one-body limits. It is 
readily seen that the original kernel K may be eliminated in favor of the reduced one, U, 
so as to obtain the following 4D equation for the t-matrix, which is equivalent to Eqs. ( |1 . 1|) 
and (P|), 

T = U + UQ T, (1.9) 

where 



g = G + G c . (1.10) 

We are now in a position to remedy the defects of our previous three-dimensional reduction. 
Applying the same 3D reduction procedure as above to this new 4D equation produces 

T 1 = U 1 + U 1 {G )T 1 , (1.11) 

where the 3D propagator is 

i Q (P)) _ Aj-Aj- A^As ATAJ A{A^ 

\ y ^)> po_ ei _ £2 24-P° + e 1 + e 2 p0-24 + e 1 + e 2 P° + e 1 + e 2 ' 1 } 

and ^2 is a parameter that arises from the eikonal approximation. This 3D propagator is 
found by integrating the four-dimensional propagator over the time component of relative 
momentum. In configuration space this action is equivalent to considering the propagator 
in which the two particles involved are considered only on equal time slices. Thus this will 
be referred to as the ET propagator. It was derived for use with instant interactions by 
Mandelzweig and Wallace 0,111 with the choice k§ = P°/2 - (ml - m|)/(2P°). It has the 



correct one-body limits as either particle's mass tends to infinity and has an invertible form. 
Meanwhile, the 3D kernel U\ is systematically defined by two equations that are obtained 
from Eqs. (|1.4j) and ( |1.5|) by the substitutions K — > U, K\ — > U\ and Go — ► Go- 

Equation ( |1 . 1 1| ) is a 3D scattering equation that incorporates relativistic effects system- 
atically and has the correct one-body limit. Numerical calculations by Nieuwenhuis and 
Tjon [||], and in our previous paper suggest that the three-dimensional equation with a 
lowest-order kernel provides a good approximation to the relevant physics of the full scat- 
tering series of ladders and crossed ladders. 

Negative-energy states and thus Z-graphs, such as the one shown in Fig. [I], are included 
in the 3D formalism in a way that is symmetrical with respect to interchange of particle 
labels. When the terms in the propagator involving Ag are omitted, it takes a form that is 



very similar to the 3D propagator of the spectator formalism of Gross [12|, as derived with 
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particle 1 on mass shell. Correspondingly, when the terms in the propagator involving A 1 
are omitted, the propagator is very similar to that of the spectator formalism, derived with 



particle 2 on mass shell. Indeed, in ++ sectors, Eq. ( 1.1 2|) is the same as either of these 
spectator propagators, and thus differences arise only because of the negative-energy states. 
Usually for the NN system, a symmetrized form of the spectator propagator is used, and 
this is obtained by averaging the propagators derived with particle 1 and particle 2 on mass 
shell [I3|. However, because of averaging, the symmetrized spectator propagator has half the 



propagation amplitudes of the ET propagator in H — and — h states. A similar comparison 
may be made with the BSLT quasipotential propagator of Blankenbecler- Sugar |I4J and 



Logunov-Tavkhelidze |7] that has been used in the work of Hummel and Tjon |15HI7[. The 
BSLT propagator also has half the propagation amplitudes of the ET propagator in H — and 
— h states. Because similar couplings to the H — and — h states are present in all of these 
approaches when similar meson exchange interactions are used, we would expect the role of 
negative-energy states to be larger when the ET propagator is used that when either the 
spectator propagator of Gross or the BSLT propagator of Hummel and Tjon is used. 

A straightforward analysis of time-ordered perturbation theory graphs for the AW t- 
matrix in the static limit (see Appendix |A|) shows that the leading Z-graphs involving the 
intermediate H — and — h states are correctly given by the ET propagator. Results based 
upon either the symmetrized Gross propagator or the BSLT propagator are too small by a 
factor of two. 

Our preference for the ET formalism is based on three facts: it embeds the correct one- 
body limit for either particle as the mass of the other particle tends to infinity; it provides the 
correct isoscalar Z-graph contributions to leading order in 1/M for the AW system; and the 
systematic 3D formalism associated with it includes retardation effects without possessing 
the unphysical singularities that generally are present in quasipotential theories. 

For relativistic bound states, the 4D vertex functions at total momentum P obey 

r(P) = U(P)Go(P)T(P). (1.13) 

Equivalently, one may use the 3D vertex function that is obtained with the use of the ET 
propagator, 

r ET (p) = u 1 (p)(g (p))r ET (p). (i.i4) 

Shortly, we shall need the relations of 4D vertex functions to the 3D ones, which is as follows, 

r(p) = q r (p)t et (p), (i.i5) 

where Q R is defined by, 

n R (p) = i + \u{p) - u x {p)\g Q {P)n R {p). (i.i6) 

A similar relation is needed for final states: 

f(P) = f ET (P)fi L (P), (1.17) 

n L (p) = i + n L (p)g (p)[u(P) - u^p)}. (lis) 
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In the spirit underlying the use of an effective NN interaction, the 3D interaction XJ\ is 
truncated at "lowest order" within a systematic expansion in powers of the coupling con- 
stant. This truncation violates Lorentz invariance, although the full 3D reduction formalism 
is equivalent to the 4D formalism and thus respects Lorentz invariance. For electromag- 
netic matrix elements, the absorption of the virtual photon's momentum Q causes the final 
state to have nonzero three- momentum, even if the initial state did not. Thus, in general, 
a dynamical boost is needed to obtain wave functions in frames where the bound state has 
nonzero three-momentum. In our 3D formalism, the dynamical boost is embedded within the 
equations. Wave functions corresponding to nonzero three-momentum should be obtained 
straightforwardly by solving the bound-state equation with the interaction appropriate to 
the moving frame. This differs from the interaction in the cm. frame because retardation 
corrections and Dirac spinors in the matrix elements depend upon the total momentum. 
The violation of Lorentz invariance causes the total bound-state energy -E'(P) to differ from 



+ P 2 . For the deuteron, the effect is small and it may be compensated by a simple 
renormalization of the interaction that, in effect, approximately takes into account the terms 
omitted when \J\ is truncated at lowest order. Results for deuteron form factors have been 
found in previous calculations to be essentially the same whether the interaction is renor- 
malized or not, thus indicating negligible sensitivity to the truncation and the associated 
violation of Lorentz invariance. 

In order to confront the predictions of this formalism with electron-scattering data we 
must derive a 3D reduction of the electromagnetic current that is consistent with the reduc- 
tion of the scattering equation. A clear 4D formalism for the current follows from coupling 
photons everywhere in the Bethe-Salpeter Green's function. This leads to the following 
gauge-invariant 4D result for the photon's interaction with the two-nucleon system: 

A BS ,» = T(P')G2 ifl T(P) + T(P')G (P')K;(Q)Go(P)T(P), (1.19) 

where P' = P + Q, (jq,« is a five-point function for insertion of a photon of momentum 
Q and Lorentz index /i in the free propagator G , and K 1 is the five-point, two-nucleon 
irreducible NN — > NN'y amplitude found by coupling the photon to all charged-particle 
lines inside the Bethe-Salpeter kernel K. A similar result holds for the 4D theory involving 
Gc-, 

Am* = r(p')g^r(p) + f (p')g (p')u;(Q)g (p)T(p), (1.20) 

where Qq^ = Gq ^ + G C41 , and UJ is the result of coupling the photon to all internal lines of 
the reduced interaction U. In order to develop a conserved current in the 4D theory based 
on Eq. ( |1.9|) , we must construct G C411 which is the extra current required in order that Qq ^ 
satisfies a Ward-Takahashi identity involving Q . Unfortunately, this is not straightforward, 
since in general the variable n\ in Gc depends on the three-momentum of the intermediate 



state. As we shall see in Section [TV], we cannot simply use the standard fermion electromag- 
netic current, 7 M , since this choice violates current conservation. In Section [TV] we construct 
the additional pieces of the current that restore current conservation, thereby obtaining a 
conserved current Qq„ for use with vertex functions obtained from Eq. ( |1.13| ). 

In this paper, we develop an equivalent 3D form of the current matrix elements for the 
equal-time formalism, based on an expansion of the formula, 
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a et ,, = f ET (p')(^ L ( J p / )^^(^)) r ET(P) + r ET (p')(g(p')u;g(p))T ET (p). (1.21) 

The resultant reduction is consistent with the reduction used to obtain the 3D t-matrix of 
Eq. Q1.11D and the vertex function of Eq. (|1.14j ). Furthermore, this reduction preserves the 
two-body electromagnetic Ward-Takahashi identities which are present in the 4D theory. 
Quantities inside angle brackets are integrated over time-components of relative momenta, 
thus reducing them to a 3D form. This yields a consistent three-dimensional formalism 
that includes the effects of relativity systematically, has the correct one-body limits, and 
maintains current conservation. We then apply this machinery to the calculation of electron- 
deuteron scattering in the impulse approximation. 



Previous calculations of electron-deuteron scattering by Hummel and Tjon |15Hl7|| have 
used instant interactions and a form of the ET current. However, several approximations 
were employed, such as the use of wave functions based on the 3D quasipotential propagator 
of Blankenbecler-Sugar [|14[] and Logunov-Tavkhelidze J7|, approximate boost operators, and 
an electromagnetic current which was conserved only in positive-energy states. Calculations 
of elastic electron-deuteron scattering also were performed by Devine and Wallace using an 
instant interaction and the ET propagator of Mandelzweig and Wallace, Eq. (|1.12|) with 
the choice of given below that equation. In that work, a suitable conserved current 



was derived for use with vertex functions that included negative-energy components. In this 
paper, we extend these previous analyses by use of our systematic 3D formalism with the 
improved choice of k\ that is needed when retardation effects are included. 

It is easy to show that if the instant approximation is used to reduce the four-dimensional 
equation, Eq. ( |1.9| ), to the three-dimensional ET equation, a particularly simple three- 
dimensional conserved current exists. If we denote this current by G\ nst ^, then we may 
write the resulting Ward-Takahashi identity as: 

Wnst>i,p 2 ; Q) = <&[<f?o>(pi, P2; P°) - (SqXpi, p 2 + Q; P° + Q )} 

+ (l<->2), (1.22) 

where (1 2) indicates the part of the current proportional to charge qi, with momenta and 
labels of the two particles interchanged. Consequently, if we construct the current matrix 
element 

Anst, M = / inst ( P ; P + Q)GLtA^ P ' Q) r inst(p; P), (1-23) 

where r inst is the vertex function obtained from the bound-state equation corresponding to 
Eq. (|1.14| ) with an instant interaction chosen for Ui, then the result will be a conserved 
impulse current. A similar instant analysis was performed by Devine and Wallace [ITS 



Our instant calculations differ from those of Ref. |19| in that a slightly different choice for 



is made than was made by Wallace and Mandelzweig. This different choice for k° is 
necessitated by the requirement that the three-dimensional theory with retardations should 
not possess any unphysical singularities when evaluated in a frame where the two-body 
system has large momentum. 

A first step beyond the instant approximation may be obtained by replacing the vertex 
r inst , which is calculated with instant OBE interactions, with the vertex Tet, calculated 
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with the full retarded OBE interaction obtained in the systematic ET formalism. In or- 
der to maintain current conservation our approach demands that in making this change we 
should also replace Qi m t,u by a significantly more complicated object, Qq^. The resultant 
current matrix element, which we denote *4et>> differs from Eq. (|1.23|) by replacement of 
the subscript inst by the subscript ET in the vertex functions and replacement of ^in St)M 
by Qo t u- This *4et> would satisfy current conservation when retardations are included pro- 
vided that we also included the two-body currents, such as the one depicted in Figure 
that become necessary because of retardation effects. This work is concerned primarily with 
displaying the formalism and performing impulse approximation calculations that explore 
the effects of retardation and negative-energy components in the vertex functions. There- 
fore, we use the simpler current (Qi nst4l ) with the ET vertex functions as well as with the 
instant vertex functions. For elastic electron-deuteron scattering this is expected to be a 
very good approximation to the use of the full ET current. Calculations where (G c m ) 
is omitted from the instant current matrix element suggest that the total contribution of 
(G G ) to the observables is small. This is consistent with the interpretation of (G G ) as 
arising from the coupling of the photon to internal charged lines in the crossed-box graph: 
(Gq ) is a two-body current, and therefore we would expect it a priori to be small at the 
momentum-transfers under consideration here. Since (G Cll ) itself makes a small contribu- 
tion to observables, using a form of it that only approximately satisfies current conservation 
is expected to have an even smaller effect on our numerical results. We expect that these 
effects, and the other effects of retardations we have neglected in the current, will only have 
a minor influence on observables. Here our main goal is to examine the effect of retardation 
on the vertex functions, and the role that negative-energy states play in the calculation, and 
so we do not calculate these additional effects. Work is in progress to include these extra 
terms of the current that are derived in Sec. [TV] in the calculation. We are also calculating 
meson-exchange currents, such as p7c~f and uja'j, which we do expect to influence observables. 

The first step in obtaining theoretical predictions for the experimental observables in 
electron scattering on the deuteron is to find the vertex functions r inst and Tet- We take 
the four-dimensional kernel U in Eq. ( |1.9|) to be the sum of single-boson exchanges. The 
parameters for these exchanges are taken from the Bonn-B OBE model, with the exception 
of the a-meson coupling, which we leave as a free parameter. From the 4D kernel, t/^oBE; an 
instant interaction, L^ingt, is easily found, and a corresponding retarded interaction, U^et, 
can be defined by the systematic procedure outlined above. These two different interactions 
are then inserted into the bound-state equation derived from Eq. ( |1.14| ). The u-meson 
coupling is adjusted so that the NN bound-state pole in the 3 Si — 3 D 1 channel appears at 
the deuteron mass. Once this is done the deuteron vertex functions r inst and T ET can be 
extracted. We calculate these two vertex functions including the effects of negative-energy 
states and also in the approximation where only positive-energy states contribute. 

With vertex functions in hand and using the 3D current (Qi nst ^), we calculate the 
electron-deuteron scattering observables A, B, and T 20 . The results from such a calcu- 
lation are shown in Figs. ^ and We also show experimental data from Refs. 21 
for A, from Refs. [^.^3|.^,|(| for B and from Refs. |27H32| for T 20 . It should be pointed 
out though, that a number of two-body effects, such as the differences between Gi nst ^ and 
Qq^ that are needed to restore current conservation, and the usual pttj MEC contribution, 
should be added to our calculations before they can be reliably compared to experimental 
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data. One of the most interesting features is the close similarity of the results based on 
the ET vertex functions that include retardations and instant vertex functions that do not. 
We find that once the a coupling is renormalized at Q 2 = (as must be done to refit the 
deuteron binding after incorporating the repulsive effects arising from meson retardation), 
the deuteron properties in these two models are remarkably similar. 

For comparison, in Figs. || |7] and [8] we display calculations for A, B and T20, where the 
effects of negative-energy states are removed. We see that the inclusion of these states in 
the calculation makes little difference to any of the observables. However, comparing the 
different positions of the minima in Figs. |] and |7] we see that including the negative-energy 
states in the calculation does shift the minimum in B to somewhat larger Q 2 . A similar 
effect was observed by van Orden et al. [3^] in calculations of electron-deuteron scattering 
using the spectator formalism. However, note that here, in contradistinction to the results of 



Ref . |33| , the inclusion of negative-energy states does not bring the impulse approximation 
calculation into agreement with the data. 

The fact that negative-energy states seem to have a smaller effect on observables in the 
ET analysis than in the spectator analysis of van Orden et al. |!3| is somewhat surprising. 
As pointed out above, the ET propagator has twice the negative-energy state propagation 
amplitude of the Gross propagator. Thus, other differences between the ET and specta- 
tor models, not just differences in the role of negative-energy states in the two approaches, 
appear to be responsible for the differing results for B(Q 2 ). The fact that Ref. [f^jfl's cal- 
culation essentially agrees with the experimental data for B(Q 2 ) is not solely attributable 
to the Gross formalism's treatment of negative-energy states, since our results show that 
negative-energy states have a much smaller effect. 

We also note that the inclusion of retardation moves the minimum in B a little higher 
in Q 2 but otherwise has little observable effect. Again, this additional effect is not enough 
to bring the predictions of our model into line with the experimental data. 

Finally, examining the tensor polarization T20 we see that all the different models produce 
results which are very similar. This suggests that this observable is fairly insensitive to 
dynamical details of the deuteron model, at least up to Q 2 = 4 GeV 2 . 

We find that elimination of the approximations used by Hummel and Tjon produces 
only minor changes for the experimentally-measured quantities A, B, and T 2 o, although the 
precise location of the minimum in the magnetic form factor B does change when meson 
retardations are included in the calculation. Nevertheless, there is little improvement in the 
agreement of the impulse-approximation calculations with the experimental data for B. 

Comparing our results with those of a non-relativistic impulse approximation calculation 
of the same observables (see, for instance, Ref. strengthens the conclusion of previous 
authors that neither the consideration of relativistic kinematics for the nucleons and mesons 
nor the inclusion of negative-energy state effects improves the agreement with the experimen- 
tal data. Of course, by definition, such a non-relativistic impulse approximation calculation 
neglects both relativistic effects and two-body current contributions. Indeed, we see here 
that, as already found by Arnold et al. [0] and Zuilhof and Tjon |5SJ, including relativistic 
effects actually worsens the agreement with the experimental data. The inclusion of these 
effects to all orders in a p/M expansion, as done here and in Refs. p4]j35| does not lead to a 
small correction which brings the theory into closer agreement with the experimental data. 
This suggests that the comparative success of a simple non-relativistic impulse approxima- 
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tion calculation is fortuitous. Dynamical mechanisms beyond the impulse approximation 
appear to play a more important role in this reaction than one would conclude from the 
non-relativistic calculation. Figures [5[ f|, and || imply that once Q 2 gets above about 0.5 
GeV 2 two-body current contributions become important. One example of such a two-body 
current would be meson-exchange current contributions, but it is also possible that other 
two-body currents are dynamically important in this regime. 

The rest of this paper is devoted to a detailed description of our formalism, and to 
the explanation of technical details pertaining to our calculation that were omitted in the 
brief sketch of the calculations given here. In Section |TJ we present our modified four- 
dimensional equation for Gc with the correct one-body limit and incorporating a form of 
the eikonal approximation which is an improvement over that proposed in Ref. [jlj. The 
new choice of the parameter k° has been found to be necessary, when retardations are 
included in the ET interaction, in order to avoid unphysical singularities when P is large. 
Our formalism for reductions to three dimensions is then used to determine the lowest- 
order 3D interaction U\. In Section |TXT| we explain the various interactions that are used 
in calculations of deuteron wave functions. These can be divided into two classes: instant 
interactions, and ET interactions that include meson retardation. Within either of these 
classes, versions of the interactions are constructed that do and do not include the effects of 
negative-energy states, in order to display the role played by such components of the deuteron 
wave function. Section |TV| discusses our 3D reduction of the electromagnetic current that 
maintains current conservation, including a detailed discussion of the way in which the 
gauging of the propagator Gc in our modified four- dimensional equation is performed. The 
ET current Qq ^ is derived and we discuss its connection to the simpler current G^st u that 
is used in the present calculations. Section presents details concerning the application 
of this formalism to the calculation of electromagnetic observables using the wave functions 



computed as discussed in Section [II. Conclusions are presented in Section VI 



II. BOUND-STATE EQUATIONS WITH CORRECT ONE-BODY LIMITS 

A. A four-dimensional equation 

As outlined in the Introduction, a simple four- dimensional equation may be obtained by 
employing a form of the eikonal approximation to reorganize the generalized ladder Bethe- 
Salpeter kernel, based upon the fact that we may always write the kernel K in the iterative 
form, Eq. ( |1.8| ). The question is: what are good choices for Gc and U? To answer this 
question, consider the lowest order (in coupling constant) parts of the kernels K and U. 
Obviously, the second-order pieces must be the same, i.e., 

U^=K^. (2.1) 

Meanwhile, at fourth order, we obtain 

=UW+U<®GcU®. (2.2) 

Thus, in order that the expansion for U can be truncated at second order without losing 
the one-body limit as m<i — > oo, the choice of Gc must be such that the last term captures 
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the part of that is non-zero in that limit. Moreover, if Gc is chosen in this way, the 
one-body limit will then be correct no matter the order at which the expansion of U is 
truncated. 

Consider the fourth-order contribution to K, i.e., the crossed-box graph depicted (with 
momentum labels) in Fig. |9|. In order to express in the form (pD, it is necessary 

to commute vertices until they appear in the same order as in K^GcKW. In doing this 
commutator terms are collected in U^. Once this is done, the part of the graph that takes 
on an iterative form is expressed as 

K^ I (k' 1 ,k' 2 ;k 1 ,k 2 )=t J ^K^(k' 2 -p 2 )d l (P-p 2 )d 2 (k 2 + k' 2 -p 2 )K^(p 2 -k 2 ). (2.3) 

Here K^> is the ladder Bethe-Salpeter kernel, P is the conserved, total four-momentum: 

P = k 1 + k 2 = k[ + k' 2 , (2.4) 
and the propagators d{ are given by, 

di(Pi) = — o 7 r— ; : n ' 7 s — > (2-5) 

Pi ~ e i(Pi) + ir l Pi+ e i(Pi) ~ l V 



where 



e t ( Pl ) = (pl + mi)V\ (2.6) 

for spin-zero particles, 



Af(Pi) { ^fezm±!^ for spin-haJf particles, ^ (2 '° 

and 7] is a positive infinitesimal. 

In Eq. ( |2.3| ) the argument of the function d 2 may be rewritten: 

(*° + A#-p°,k 2 + k4-p 2 ). (2.8) 

Suppose that the zeroth component of the momentum of particle 2 is large. In particular, 
this will be the case if m 2 ^> mi (one-body limit). In this case particle two's intermediate 
and final-state three-momentum will be largely unaffected by the presence of particle one, 
and so we may approximate these components as unchanging, 

k 2 + k' 2 w 2p 2 . (2.9) 

This is the eikonal approximation, and it rests upon forward scattering being the dominant 
mechanism of interaction. Indeed, making the replacement ( |2.9| ) in ( |2.3| ) will not affect the 
value of i^tr m the limit m 2 —>■ oo. 

This argument shows that K^ r may be approximately rewritten as: 

K^(k' 2 ,k 2 ;P)Kij ^K^(k' 2 -p 2 )d 1 (P-p 2 )d 2 (2n 2 -p 2 )K^(p 2 -k 2 ), (2.10) 
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with the four-component object k 2 defined by: 



In operator notation 

4tl ~ ^ (2) G C ^ (2) . (2.12) 

Now we are in a position to answer the question of how to choose Gc- Equation (|2.10|) 
shows that the choice: 

G G (p;P)=id 1 {is 1 P+p)d 2 {2K 2 -is 2 P+ P ), (2.13) 

which we write symbolically as, 

G c = id 1 d c 2 , (2.14) 

will allow U to be truncated at second order while still yielding a four-dimensional equation 
with the correct one-body limit. Here it is understood that when writing Gq as a function 
of the four four- momenta p[, p' 2 , pi, and p 2 , we have 

Gc(p' 1 ,P 2 ;Pi,P2) = i(2n) 8 S^(p[-p 1 )S^(p' 2 -p 2 )d 1 (p 1 )d 2 (2 K2 -p 2 ). (2.15) 

On the other hand, in writing Eq. ( |2.13| ) we have expressed the result in terms of center-of- 
mass and relative four-momenta: 

p 1 = v 1 P + p, Pi = v 2 P-p\ (2.16) 

P 2 +m\ — ml P 2 — m\ + ml 
Ul = 2P 2 ' " 2 = 2P 2 ' ^ 2 ' 11 > 



The propagator d 2 defined in (|2.13|) depends on the zeroth components of the external mo- 
menta through the zeroth component of k 2 . Thus the use of operator notation in Eq. ( |2.12p 
anticipates a suitable choice for these zeroth components of external momenta. 

In the case that particle two is extremely massive, it will not be very far off its mass 
shell. The total energy P° of the system will be approximately the energy of particle two, 
and thus we choose 

nl = u 2 P°. (2.18) 

This choice was used in Ref. |j| because it leads to an integral equation that satisfies both 
the criteria needed to obtain the correct one-body limit m 2 —>■ oo. It also is suitable for 
equal-mass particles in the center-of-mass frame. However, in frames where the bound state 
has nonzero total momentum, the choice leads to an interaction that possesses unphysical 
singularities. 

A consideration of the singularities of the interaction when the particles have equal masses 
suggests that the particles should be equally far off their mass shells in the intermediate state. 
Hence, in that case we choose 
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-(P c 

9, V 



ei + e 2 ), (2.19) 

where, in an arbitrary frame, 

ei=ei(P/2 + p); e 2 = e 2 (P/2 - p). (2.20) 

This choice is consistent with the eikonal approximation and differs from that of Eq. (|2.18|) 
only when P / 0. Note that our new choice for k\ depends not only on the external 
variables P° and P, but also on the internal three-momentum p. This feature complicates 
the construction of a conserved current. 

Once U and the Gc of Eq. ( |2.13| ) are chosen, this defines an "improved" ladder BSE, 
which, in the two-body center-of-mass frame takes the form, 

d 4 p 



r(p'°,p'; S ) = i J -^-U(p',p)d 1 (u 1 P +p°,p) 



d 2 (u 2 P° - p°, -p) + d 2 (24 - u 2 P° + p°, -p)] r(p°, p; s). (2.21) 

Because the choice of employed here agrees with that of Ref. in the limit m 2 — > oo the 
proof given there now suffices to show that this equation becomes the appropriate one-body 
equation in the limit m 2 — > oo. For equal mass particles, is chosen as per Eq. ( |2.19|) . 
We make the choice ( |2.19|) because it eliminates singularities that afflict the choice Q2.18|) . 



As observed in Ref. jl] if the choice ( |2.18| ) is made for k\ then the equation ( |2.21| ) develops 
a cut beginning at 

2z/ 2 P° = 2m 2 + ii. (2.22) 



If we attempt to solve the integral equation (|2.21| ) in a frame moving with large enough three 



momentum P, then P° = JMj + P 2 will become big enough to encounter this singularity. 
However, with the choice ( |2.19| ) the interaction does not have this unphysical singularity. 

The above arguments pertain to the one-body limit in which m 2 — > oo. Clearly we can 
interchange the labels on particles one and two in order to obtain a propagator id\d 2 that 
is motivated by a consideration of the limit m\ — > oo. In the propagator d\ we choose, for 
the equal-mass case 

K o = I(P° + ei-e 2 ). (2.23) 



In the case m\ = m 2 , writing the propagator as in Eq. ( |2.14| ) treats the two particles 
asymmetrically. If they are identical particles this will lead to violations of the identical- 
particle statistics. To avoid this, we symmetrize the propagator Gc, by choosing, instead of 
idid 2 , the form 

Go = Gq + Gc 

= t -(d 1 + d c 1 )(d 2 + d c 2 ). (2.24) 

If the mass of particle n is taken to infinity we use the choice k° = P° in d c n ^ while the k° 
variable associated with the light particle is chosen to be zero. This yields the one-body 
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propagator for the light particle, multiplied by a delta function in the relative momentum, 
and a projection operator onto the positive-energy states of particle n. Meanwhile, the delta 
function guarantees that the interaction will be static, and so both conditions necessary 
for the correct one-body limit to be present are satisfied. However, note that this formal 
property does not apply to the case of practical interest for this work, mi = m 2 , since then 
k§ and k° 2 are defined by Eqs. (ggg) and ( pTTSf ). 

Symbolically we now write our "improved" ladder BSE as 

r = ug r, (2.25) 

where here, and throughout this work, Qq is defined by Eq. ( |2.24 ). 

We stress that what has been done here is to take certain pieces of the Bethe-Salpeter ker- 
nel K and rewrite them as K^G C K^ 2 \ K^G C K^G C K^ 2 \ etc. Consequently, Eq. (^25|) 
is equivalent to a Bethe-Salpeter equation in which graphs other than one-meson exchange 
are approximately included in the kernel. Thus we expect that the solution of this equation 
may provide a better description of the dynamics of two-particle systems than the ladder 
BSE amplitude. Numerical calculations in a scalar field theory appear to support this [p]. 



B. Summary of the reduction to three dimensions 

The 3D reduction that was outlined in the Introduction may now be applied to the 4D 



equation ( 2.25 ). This yields the following 3D equation for the bound-state vertex function, 
Tet, 

r E T = tA(£o>r E T, (2.26) 



with (Q ) the three-dimensional propagator of Eq. ( |1.12|) , and the interaction U\ defined as 
the solution of the following coupled equations, 

Ui = (g )- 1 (g ug)(g )-\ (2.27) 

which is three-dimensional, and 

g = g + g (u - u x )g. (2.28) 

Note once more that the solution for U\ is found by demanding that (g) = (go)- Equa- 
tions (p.26|) and (|2.27|) are exactly equivalent to the 4D equation ( |2.25| ). Inserting the choice 



( p,19| ) for k\ and specializing to the case of equal-mass particles gives 

(c (P)] A ^ + A ^ A 2 A rA 2 + ArA 2 - 

which is a restatement of Eq. (|1.12j) . We note in passing that this is exactly the same 
expression as that found when the simpler expression i(c?i(c? 2 + <#D) * s calculated. 

By using Eq. (|2.26| ), rather then Eq. (|1.2|) , we not only have a consistent 3D reduction 
but we approximately include the crossed-ladder graphs, particularly the "Z-graph" contri- 
butions where one nucleon is in a positive-energy state and the other is in a negative-energy 
state. 
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For equal- mass particles in the center-of-mass frame, Eq. (|2.26|) agrees with the quasipo- 
tential equations of Wallace and Mandelzweig [11|. However, in a quasipotential approach 
the interaction in other frames should be obtained by boosting via a dynamical equation. 
This leads to unphysical singularities. The systematic ET formalism avoids these. In essence, 
the ET formalism with a truncated interaction emphasizes the elimination of unphysical sin- 
gularities over the strict enforcement of Lorentz invariance. 

(2) 

The second-order 3D interaction U{ is given by: 

(Qo)ui 2) (g ) = (GoK^g )=A. (2.30) 
An explicit form for A can be computed for the case of OBE interactions of the form 

X»(,)=lM (2.31) 

n Q Pn 

where V n (i) denotes the appropriate vertex operator for the interaction of the nth meson 
with nucleon i. If we write: 

\P1\P2 API AP2 

A = Y1 E -^AMP^ PiP2)V n (l)Vn(2)^-^, (2.32) 
n PlP > lP2P > e i + e 2 ei + e 2 

with ei = piViK® — £j, e\ = p'^i^i — 4; then, for the exchange of meson n, the factors 
A n (p'iP2 *~ P1P2) can be written in any rho-spin channel as 

A n (p', p; P°, P) {p' lP ' 2 <- p lP2 ) = ^- \F(u n + q°) + FK- q°j\ , (2.33) 



where 



r,/ X 1 1 1 1 



ei + e 2 — uo e'i + e 2 — 10 t\ + e[ — uo e 2 + e' 2 — uo 
_^ I / ei + e 2 V 1 I ei + e 2 \ 



+ 



ei + e' 2 — to \ e 2 + e' 2 — uo J e[ + e 2 — 00 \ e\ + e[ — uo / 
/ ei + e' 2 \ 1 / e\ + e' 2 \ 1 



ei + e[ — uo J e\ + e 2 — uo \e 2 + e 2 — uo J e[ + e 2 — uo 



+ e '^ 4 ( ' + ^ ) ei+ ; 2 , (2.34) 

e 2 + e 2 — u; \ei + e 2 — e x + e 2 — u; / e\ + ei — c<j 

and go — — K i- The key difference here from the calculation of Ref. |J is that here 
the parameter k® varies with the intermediate-state momentum under consideration. For 
equal-mass particles 

" L V + ei -e 2 ), 4 = l(P°-e 1 + e 2 ); (2.35) 



1 2 v ■ -x - v , z 2 
<° = J (P + ^ - 4) , 4° = 4 + 4) • (2-36) 
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The expression for A is easily converted into a result for Ui using the inverse of the 
propagator (Go). The result is expressed in terms of matrix elements between initial and 
final states of any p-spin as follows, 



<^' p1P2 (p',p; p°,p) = E [^ViPiKaKHpiPi)] [^(p 2 p 2 )v;(2)^(p 2 p 2 ) 

xA n (p',p; P°,P)(p , 1 p 2 ^p 1 p 2 ) 



(2.37) 



This is the full ET interaction including retardation effects. If one eliminates the retardation 
effects, the result is the instant interaction defined by, 



Y, [^ViP'iKaKKpiPi)] [^V 2 p 2 )K(2K 2 (p 2 p 2 



-9n 



Expanding the vertex functions in terms of Dirac spinors, 

r ET ( P ,p) = 7 ?7 2 ° E « pl (piPi)« P2 (p2P 2 )r^ 2 ( P ,p), 

P1P2 

leads to the coupled equations that we solve, 

^ d s p' 



(2.38) 
(2.39) 



- IP-P) = E / ^^^(P,P';^,P)<ft)> pi/,& (p',P)I^(p',P)- (2-40) 



-1P1,P2/ 



Here the propagator factors (So)"' 1 " 2 are (£ > ++ = l/(P°-e' 1 -e' 2 ), (£ > + ~ = l/2e' 2 , (£ >~ 



l/2ei and (£ >~ 



-l/(P° + e ' 1 + e 2 ; 



Two special cases of particular interest are the matrix element of the ET interaction 
between positive-energy spinors, as follows, 

C/S ++ (p',p;P Q ,P) = E [^(p'i)K(lK(pi)] [u + {v' 2 )V n {2)u + { V2 



X- 



2coV, 



1 

+ 



P°-e 1 



e 2 



po _ e / _ e2 _ ^ n 2 (P° - e' a - e 2 - ^ n ) 2 

, 1 ^-el-e, l (P°-6;-6 2 )(P - ei - e2 ) 
2 (P° - ei - e 2 - cu n ) 2 2 (P° ' ! 



e' x - e 2 - ^n) 3 



f2.41^ 



and the matrix element in positive-energy states omitting the effects of Gc, 
^6pt(^,P;^°,P)= E [« + (p!)K(1)« + (Pi)] [n + (p 2 )K(2)n + (p 2 ) 

' +(1-2 



A 

2uj r , 



x 



po _ e / _ e2 _ ^ 



^2.42) 



This is the standard time-ordered perturbation theory (TOPT) one-boson-exchange inter- 
action. Diagrams for the one-pion exchange part of this interaction between positive-energy 
nucleon states are shown in Fig. [10]. There are additional pieces in the ET interaction of 
Eq. Q2.41 ) that arise from Gc, i-e., from the approximate inclusion of higher-order graphs 
in the crossed-ladder kernel, which is necessary to obtain the one-body limits. 
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III. CONSTRUCTION OF THREE-DIMENSIONAL NUCLEON-NUCLEON 

INTERACTIONS 



A number of different 3D two-body interactions may be defined by including or not 
including retardation effects, and by including, or not including the contributions of negative- 
energy states. These different 3D interactions are used in our assessment of the significance 
of retardations and of negative-energy states in calculations of electromagnetic observables. 
In this section we outline the development of the various interactions used in this work, 
in each case starting from a four- dimensional kernel K that is the sum of six single-meson 
exchanges. The mesons are the tt(138), the <r(550), the ^(549), the p(769), the u;(782), and 
the 5(983). The quantum numbers and masses of these mesons, the cutoffs in the (monopole) 
form factors, as well as the couplings for all but the a meson, are listed in Table |. All these 
parameters except for the a coupling are taken directly from the Bonn-B fit to the NN 
phase shifts p6| . The o coupling is varied so as to achieve the correct deuteron binding 



energy for each interaction considered. 

We then use the techniques of the previous section to construct the following interactions, 
all of which are to be used in the three-dimensional ET equation (|2.4U|). (Note that in 



calculating the 3D interaction we assume that the 4D interaction contains no form factors 
and any dependence of vertex factors V n upon q° is neglected. After calculating the three- 
dimensional interaction we then insert monopole form factors at all the vertices.) 



1. The "Retarded ET" interaction, given by Eq. ( |2.37|) , with the choices ( [2.23|) and (|2.19|) 
for Ki and k\. 

2. The instant NN interaction, defined by Eq. (|2.38j ). 

3. The second-order TOPT NN interaction defined by Eq. (|2.42|) . 



The first and second interactions are used in a two-body equation with the full ET Green's 
function given by Eq. ( |2.40| ), and also in an equation in which only the ++ sector is retained. 



However, for the instant interaction, we follow the practice of Devine and Wallace [18[] and 
switch off by fiat couplings that involve flip of both p-spins, i.e., between the ++ and — 
sectors, and between the H — and — h sectors. A partial justification of this rule follows 
from an analysis of the static limit of our 3D retarded interaction, which shows that all 

couplings to the states vanish as 1/M — > 0, while other couplings approach the instant 

form. Although the static limit for coupling of H — and — h sectors approaches the standard 
instant form, we omit this coupling in order to reproduce the previous results of Devine and 
Wallace. 

Since, as we have already discussed, the propagator (Gq) _1 does not exist in the H — and 
— h sectors, the TOPT interaction, U^topti is used with only ++ sectors retained in the 
equation. 



Once a particular interaction is chosen, the integral equation ( |2.40| ) is solved for the 



bound-state energy. The method used involves seeking the energy at which the largest 
eigenvalue of the kernel U\(Qq) is one. The eigenvalue is calculated using the Malfliet-Tjon 
iteration procedure. Details of this method, the rho-spin basis chosen, and the way the 



angular integrations are performed may be found in Ref. [19|. In our calculations we used 



40 quadratures in the radial momentum, and 8 in the integration over the polar angle, 9. In 



17 



the case of the two energy- dependent interactions numerical integration was used to perform 
the integration over the azimuthal angle, <fi. We found that 12 quadratures were sufficient 
to achieve convergence. 

In each calculation, the a coupling was adjusted to get the correct deuteron binding 
energy, producing the results (accurate to three significant figures) given in Table [H]. The 
value given for the instant calculation with positive-energy states alone is that found in the 
original Bonn-B fit. In all other cases the a coupling must be adjusted to compensate for 
the inclusion of retardation, the effects of negative-energy states, etc. We believe that this 
adjustment of the scalar coupling strength is sufficient to get a reasonable deuteron wave 
function. However, one direction for future work is to refit the meson-exchange parameters 
in these various different relativistic NN interaction models to the NN scattering data. 
From Table [Tl] we see that (with the exception of the TOPT interaction), adjusting the 
a coupling to reproduce the observed deuteron binding leads to at most a 5% deviation 
from the Bonn-B value, thereby suggesting that other meson parameters would change only 
slightly if a more detailed fit were performed. 

Once the bound-state wave function in the center-of-mass frame has been determined 
in this fashion, it is a simple matter to solve the integral equation ( |2.26|) in any other 



frame. As we shall explain below, we choose to calculate electron-deuteron scattering in 
the Breit frame. Hence, we need to calculate the deuteron wave functions in frames with 
total four-momentum (yj + P 2 , P). To do this, the interaction is recalculated in the 
new frame using the rules given above, and then the integral equation is solved with this 
new interaction. Because the formalism we use for reducing the four-dimensional integral 
equation to three dimensions is not Lorentz invariant if the potential U\ is truncated at any 
finite order in the coupling, we have calculated the eigenvalue A(P) defined by 

f/ 1 (p)(6?o(P))rET(P) = A(P)r ET (P), (3.i) 

for each of the five different interactions defined above. Since A(0) = 1, by construction, 



and P° = \J M]j + P 2 , is set in accord with Lorentz invariance, the deviation of A from one 
indicates the violation of Lorentz invariance in the interaction U\. The results of this test 



are shown in Fig. 11 



Three features of these results are worthy of note. Firstly, for elastic scattering on 
the deuteron at Q ~5 fm -1 in the Breit frame, the initial and final states have momenta 
Q/2 ^2.5 fm _1 and A deviates from 1 by only 2% - 3%. Thus for momentum transfers up to 
of order 1 GeV the violations of Lorentz invariance are actually quite small. Secondly, adding 
meson retardation to the formalism actually increases the violation of Lorentz invariance. 
This is a little surprising because the retardation effects, when expanded to order P 2 /M 2 , 
reproduce the Poincare boost operator of Ref . |37| H . We note that if the Poincare boost 



operator of order P 2 /M 2 were actually sufficient in this calculation, and effects of higher 
order in 1/M were truly negligible, then the inclusion of these retardation effects should 
remove most of the violation of Lorentz invariance present in the instant analysis. This does 
not happen. Finally, including negative-energy states in the formalism decreases the size of 
the violation of Lorentz invariance. 
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IV. CONSTRUCTION OF A CURRENT-CONSERVING ELECTROMAGNETIC 

INTERACTION 



In this section, we first review how a conserved electromagnetic current is developed 
for the 4D Bethe-Salpeter formalism, and then show how a corresponding 3D current is 
obtained from the 4D one. The emphasis is on satisfying the appropriate Ward-Takahashi 
identities. These reviews establish the procedure that is followed in subsequent subsections 
in order to construct the electromagnetic current for the 4D and 3D formalisms involving 
the propagator Gc- The tricky issue when Gc is present is the construction of the additional 
terms in the current to maintain current conservation. Because Gc enters as a part of two- 
body interactions associated with crossed graphs, the current associated with it is a two- 
body current. We construct the extra terms in the current required by the Ward-Takahashi 
identity, which is not straightforward because the variables an d ^2 that we choose depend 
on the three-momentum of the intermediate state. As we shall see below, this means that we 
cannot use the standard fermion electromagnetic current, 7 M , because such a choice would 
violate current conservation. We show that the Ward-Takahashi identity may be recovered 
by incorporating an additional term in the current. This yields the ET current referred to as 
Qq ^ in the Introduction. We then develop the corresponding 3D current for the formalism 
involving Gc- Because the ET current has parts attributable to two-body effects, we also 
develop a simpler current that is appropriate for use with instant interactions, ^ st „. This 
last-mentioned current is the one used in the calculations of this paper, which are based 
upon the impulse approximation. Although we do not use the full formalism developed in 
this section for calculations in this paper, the results presented here are pertinent to future 
calculations in which two-body currents will be included. 



A. Review of Ward-Takahashi identities in 4D Bethe-Salpeter formalism 

We begin our discussion by reviewing the WTI for the usual Bethe-Salpeter Green's 
function. Consider the two-body Green's function G that is the solution of the Bethe- 
Salpeter equation 

G = G + G KG, (4.1) 

where K is any BSE kernel. Define the two-body Green's function for the interaction of two 
free particles with a photon of fixed momentum Q via: 



Gi*(pi,P*,Q) = G ( Pl , P2 + Q)Htf\Q 2 K\ Pl )]G ( Pl ,p 2 ) + (1 

-,7(2 
r 0,M 



GtS (pi , P 2, Q) + (Pi,P2,Q). (4.2) 



Here, and throughout the rest of the paper, the notation (1 <-> 2) indicates that the momenta 
of the two particles must be swapped, and the labels exchanged. Therefore, the (1 <-> 2) 
pieces of any expression represent the photon coupling to whichever particle it did not couple 
in the first part of the expression. An explicit example of this rule is the particle one coupling 
term of Eq. which is Gf$(p 1} p2, Q) = G (pi + Q,p 2 )H^> 2 )j£W)]Go(pi,p 2 ) 



The free Green's function Gq ^ obeys a Ward-Takahashi identity. Now let G 1 be the 
Green's function for the interaction of one photon with the interacting two-particle system. 
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Note that in the two-nucleon system the charges % will include isospin operators, so care 
must be exercised in ordering charges and isospin-dependent interactions. We may write 
the following equation for a current-conserving G^, by allowing the photon to be inserted 
anywhere on the right-hand side of Eq. ( |4.1|) , and then rearranging the result: 

+(!<_> 2) + J gW^W^ G^, k' 2 ;p' 1 ,p' 2 )K;(p' 1 ,p' 2 ,p 1 ,p 2 ; Q)G( Pl ,p 2 ; fc x , k 2 ). (4.3) 

Here is found by coupling the photon to every internal charged line in the kernel K. 
(Note that there is an overall delta function 5^(k[ + k' 2 — k\ — k 2 — Q) on both sides of this 
equation.) Using the WTI for Gq^, Eq. Q4.3|) , and the explicit form of Eq. (|Q| ) we find a 
WTI for Gft, 

Q^Gl(k[, k' 2 ; k u k 2 ; Q) = q 2 G(k[, k' 2 -Q;k u k 2 ) - G(k[, k' 2 ; k x , k 2 + Q)q 2 + (1 <- 2), (4.4) 
provided that 

q 2 K{p' x ,p' 2 - Q;pi,p 2 ) - K(p' 1 ,p' 2 ; Pl ,p 2 + Q)q 2 + (1 <- 2) = Q^K;(p' 1 ,p' 2 ;p 1 ,p 2 ; Q), (4.5) 

which is the WTI for the interaction current. (Similar identities are used in the construc- 
tion of Ward-Takahashi identities for the Gross — or spectator — formalism in Refs. [|38| ,|39H . ) 
The result (|4.5| ) is completely general, and will always hold if the two-body current K 1 is 
constructed by coupling the photon to every charged line in the kernel K. 

Using the usual decomposition of the two-body Green's function into a pole and a regular 
part gives the amplitude for interaction of the bound-state with a photon of momentum 
Q |^(J. Expressing the result in terms of total and relative four-momenta yields: 

A,(P,Q)= J -^t{p-Q/2 ] P')Gl^{p ] P ] Q)T{p ] P). + {\^2) 

+ / ^^T(p';P')G (p';P')K;(p',P';p,P;Q)Go(p;P)T(p;P), (4.6) 

where P' = P + Q and T{p\ P) is the two-body vertex function corresponding to the bound- 
state at four- momentum P 2 = M 2 . From the WTIs for Gq^ and KJ, and the bound-state 
BSE, T = KGoT, it is then straightforward to show that 

W,Q) = 0, (4.7) 

as required for current conservation. 



B. Gauge invariance in a 3D reduction of the Bethe-Salpeter formalism 

In the context of a reduction to three dimensions, the question is how to maintain gauge 
invariance when the reduction is made. In Section | we outlined the 3D reduction method 
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for the bound-state vertex arising from the solution of the Bethe-Salpeter equation. This 
gives 

T 1 (P) = K 1 (P)(G (P))T 1 (P), (4.8) 

where the dependence of quantities on the total four momentum, P, is indicated. The 4D 
vertex function, T, and the corresponding 3D one, Ti, are related as follows, 

G (P)T(P) = g(P)T l (P). (4.9) 

Moreover, the interaction K\ is defined so as to obey, 

(G (P))K 1 (P)(G (P)) = (G Q {P)K{P)Q{P)). (4.10) 

Inserting Eq. (|4.9| ) into Eq. fl4.6|) , we find an entirely equivalent expression for the current, 



MM) = J ^'^^ U^pW^-p') 



(i) 



(2k) 

(Q 2 )d^(u 2 P - p)5^\p> -p- Q/2) - ij^iQ^d^P + p)5^(p' -p+ Q/2)) 

+KVP',p'; P,p; Q)} G(p, k; P)ri(k; P). (4.11) 



Because the vertex functions do not depend upon the time- components of relative momenta, 
integrations over k° and k° reduce this expression to a 3D one, which we abbreviate as 

a, = v x {p'){G{p') [J, + k;} ^(P))r 1 (p). (4.i2) 

Now, given a result for Ti obtained by systematic expansion of Ki, the amplitude A^ also 
can be expanded systematically in a way that maintains current conservation. First, note 
that we solve for K\ in accord with Eq. (|4.10| ). This leads to an infinite series for K\. If the 



condition ( |4.10| ) is imposed order- by-order in the expansion in K — Ki, the condition defines 
K\ to the same order. Truncation of the kernel is necessary for a practical analysis and we 
must ask if a corresponding 3D approximation for the current matrix element (|4.12j) exists 
that maintains the Ward-Takahashi identities of the theory. It turns out that the current 
matrix element ( U-1Q ) is conserved if Q(J^ + K^)Q on the right-hand side of Eq. M.lfy ) is 



expanded to a given order in the coupling constant and the kernel K\ used to define Ti is 
obtained from Eq. by truncation at the same order in the coupling constant. 

To effect this the right-hand side of Eq. ( |4. 12| ) is split into two pieces, one due to the 
one-body current J M , and one due to the two-body current K^. Suppose now that K\ has 

been truncated at lowest order, i.e., K 1 = k[ 2 \ and that K 1 = Kj( 2 \ Then, in the piece, 

we expand the C7s and retain terms up to the same order in - Kf\ A piece from the 
two-body current must be added to this. In that piece we stop this expansion of Q at zeroth 
order in — k[ 2 \ i.e., write Q = G . Thus, we define our second-order approximation to 
A„A^,by 

4 2) = f i( p ')Wri(P) 

+ UniGoiP'KK^iP') - K} 2) (P'))G^)r 1 (P) 

+ ^(P')(GU K{2) (P) ~ ^ 2) (P))G (P))r 1 (p) 

+ f 1 (P')(Go(P , )^ (2) Go(P))r 1 (P). (4.13) 
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To show that Af is gauge invariant, it must be contracted with the four- vector Q. The 
WTIs for Gq^, K^, together with the bound-state equation fl4.8[), can be used to show that 

+T 1 (P')q 1 (G (P)(K^(P) - ir{ 2) (P))Go(P))r 1 (P) + (1 <-> 2). (4.14) 
Thus, at second-order in the coupling constant, with Eq. ( |4.10|) expanded at second order 

(2) 

defining K{ , the corresponding amplitude for electromagnetic interactions of the bound 
state, as defined by Eq. ( |4.13|) , obeys 



Q»A® = 0. (4.15) 

It is straightforward to check that the same result holds if Eq. ( |4.10[ ) for K\ is truncated 
at fourth order, while the one-body and two-body current pieces of Eq. ( |4.11| ) similarly 
are expanded to fourth order. This defines a vertex function A^ which obeys Q^A^ 



0. Thus, truncating the kernel K\ defined by Eq. ( 4.10 ) and the electromagnetic vertex 
defined by Eq. (|4.12|) at consistent order in the coupling constant yields a current-conserving 
electromagnetic matrix element. A covariant extension of the formalism presented here is 
given in Ref. @. 

In fact, Aff includes contributions from diagrams where the photon couples to particles 
one and two while exchanged quanta are "in-flight" . These contributions are of two kinds. 
Firstly, if the four-dimensional kernel K is dependent on the total momentum, or if it 
involves the exchange of charged particles, then gauge invariance requires the presence of 
terms representing the coupling of the photon to internal lines in K. Secondly, even if gauge 
invariance does not require the presence of such terms in the four-dimensional formalism, 
terms arise in the three-dimensional formalism where the photon couples to particles one 
and two while an exchanged particle is "in-flight" . These must be included if our approach 
is to contain a WTI. (See Fig. |2| for a diagrammatic interpretation of one such term.) 

A special case of the above results occurs when retardation effects are omitted, i.e., the 
kernel K\ = Ki QSt , is chosen, and the bound-state equation ( |1.2|) is solved to get the vertex 
function Ti = Fi nst . Then a gauge-invariant current is obtained by implementing the instant 
approximation in the expression ( |4.13| ) in a way consistent with that in which it was used 
in obtaining the Salpeter equation ( |1.2| ). Taking Eq. ( 4.13| ) and imitating the derivation of 



the Salpeter equation by replacing K by K inst leads to 

A^P, Q) = r inst (P')(G2 M >r inst (P) + f inst (P')(G (P'))Kl st ^G (P))T inst (P), (4.16) 

where we have also replaced the meson-exchange current kernel K 1 by the instant approx- 
imation thereto. By definition this instant meson-exchange current kernel obeys the WTI 
(4~5"), but with the instant kernel Ki nst appearing on the right-hand side. To show that 



Eq. ( |4.16|) defines a gauge-invariant matrix element we contract with the four vector Q 11 and 
use the Ward-Takahashi identities for J M and K^ nstfl . This gives: 

QM inst , M (P,Q) = r inst (P')g 2 (G (P) - Go(P'))r inst (P) 

+f inst (P')(G (P / ) (q 2 K inst (P) - K mst (P')q 2 ) G (P))T mst (P) + (!<-> 2). (4.17) 
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Using the bound-state equation, r inst = Ki nst {Go)Ti nst , both at total momentum P and at 
total momentum P' then gives the desired WTI: 

Q^A inst ,,(P,Q) = 0. (4.18) 

Note that if the four- dimensional interaction K does not depend on the total momentum P 
then the simple current 

Anst,^, Q) = T inst (P')(Gl^(P; Q)>r inst (P), (4.19) 

is gauge invariant. 



C. Gauge invariance in the 4D formalism with Gc 

In this section we construct a conserved current for the 4D equation 

r = ug v. (4.20) 

In order to do this we rewrite the equation as two coupled equations, the first of which 
involves the Bethe-Salpeter kernel K, 

T = KG T, (4.21) 

and the second of which defines the reduced interaction U, 

K = U + UG C K. (4.22) 



From Eq. ( |4.6| ) we see that the gauge-invariant current for the photon coupling to the bound 
state will contain a piece corresponding to the coupling of photon lines in all possible ways 
inside K. This piece of the current is composed of two parts, as follows, 

K; = (1 + KG C )U;{1 + G C K) + KGq^K. (4.23) 

The Ward-Takahashi identities that are required to be satisfied involve one for that takes 
the same form as Eq. (|4.5|), with K replaced by U. 



q 2 U(p' 1 ,p' 2 - Q; Pl , P2 ) - U(p' 1 ,p' 2 ;p uP2 + Q)q 2 + (1 ~ 2) = Q"U;(p' 1 ,p' 2 ;p 1 ,p 2 ; Q). (4.24) 

The other involves the Green's function Gj^, representing a photon coupling to the propa- 
gator Gc- This piece must be constructed in accordance with the Ward-Takahashi identity: 

Q^GlJp^ Q) = q 2 (G c (pi,p 2 ) - G c {pi,p 2 + Q)) + (1 <- 2). (4.25) 



When both of the WTI's (|4.24j) and ( 4.25Q hold, then the electromagnetic kernel K 1 will 



satisfy the WTI ( |4.5[ ). In this case the current appropriate to a photon insertion in the free 
propagator is 

Ql^Gl^Gl^ (4.26) 
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and so, if T obeys Eq. ( f4.20|) , the amplitude 



MP, Q) = J ^-Ap - Q/ 2 > ^')< 2) (p; P QMp; p) + (i <-> 2) 



+ / ^^r(^;P0&(?/;P0^(P'.^;p.^;O)^o(p;^)r(p;P), (4.27) 

will satisfy the current conservation condition 

In order to implement this formalism, we must construct a current th at obeys 
Eq. (|4.25|) . Rather than immediately tackling the full Gc defined by Eq. Q2.24 ) we first 
perform this task for one term by considering the simple case in which the propagator is 
id\d 2 and the photon couples only to particle two. This is non-trivial because the choice of 
variable 4 in d 2 in general depends on the three-momenta of both particles. For instance, if 
4 is chosen as in Eq. fl2.19|) and the particles are of equal mass then the 4 which appears 
in the propagator Gc(pi,P2 + Q) is (for equal-mass particles) 

4' = \(Po + Qo-e 1 + e' 2 ), (4.28) 

where P = p\ + P2, and e' 2 = 62 (P2 + Q)- This is to be compared to the case where 4 is 
defined via Eq. (p. 18|) . In that case the 4 appearing in Gc(pi,P2 + Q) is, for equal mass 
particles, 

4'=\(Po + Qo), (4.29) 

and there is thus no dependence on the intermediate-state three-momenta of particles one 
and two. As we shall see below, the appearance of momenta in the choice ( f4.28| ) for k 2 leads 
to complications in the construction of a current obeying ( ]4.25| ). 
To construct such a current, we first postulate the form 

<PS ] (P2-, Q) = 4(p 2 + Q)j%d c 2 (p 2 ), (4.30) 

where we write d 2 in order to indicate that the constraint on k 2 is different once the photon 
strikes particle two. The WTI for this current is then: 

Q'jSl = Q2(4(p 2 + Q)- 1 - d c 2 (p 2 )- 1 ). (4.31) 

Constructing the right-hand side of this equation, for spin-half particles we obtain 

g 2 [2(4- fa)- Q] (4.32) 

Using the form for k 2 , Eq. ( |2.11| ), and three-momentum conservation, p' 2 = p 2 + Q, we see 
that 

4( P 2 + QY 1 - d c 2 Xp 2 y l = [2(4' - 4) - 2Q°] 7°+ Q. (4.33) 

Note that if we make the choice fl2.18j ) then, k 2 is given by 4' = ^{P + Q°) and 4 = v iP Q 1 
so 
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d\{p2 + Q)- 1 ~ dlfa)- 1 =@ + 2(v 2 - 1)Q°7 2 °. (4.34) 

This means that, in the Breit frame where Q° = 0, the usual current = q2 r )2 tl satisfies 
the Ward-Takahashi identity ( |4.31| ). 

However, when we use the formulae (|4.28 ) and ( |2.19|) in the case of equal-mass particles, 
note that 

2(4' ~ 4) = Q°+ Q " (P 1 + / P2) , (4.35) 

e 2 + e 2 

thus implying that 

4( P 2 + Q)- 1 - d c 2 { V2 )- 1 =n~Q- 4^720, (4.36) 

^2 + e 2 

where p' 2 = (e 2 , p 2 ), p 2 = (e 2 , p 2 ). Consequently, the WTI ( |4.31|) requires that we define the 
current 

jgl = ?27„-5 2) , (4-37) 



where 



= ^2— r— — 72 , (4-38) 
e 2 -+- e 2 



is the additional piece of the current that compensates for the dependence of 4 on momen- 
tum. Note that j^Q = with this choice of jj®. Note also the identity 

Oft® = 2(4' - 4)q 2l20 . (4.39) 

Next we consider a second piece of Gq, namely the propagator id\d 2l still in the situation 
where the photon couples only to particle two. Denote the free NN — > NN'-f Green's 
function corresponding to this situation by d c . The Ward-Takahashi identity which this 
Green's function should obey is, 

Q fl dl(p uP 2; Q) = iq2[d c 1 (pi)d 2 (p 2 ) - d\( Pl )d 2 (p 2 + Q)}. (4.40) 

Note that here the propagator d\(p{) ^ dl(pi), since the change in the momentum of particle 
two after the impact of the photon leads to a change in the choice of 4 that appears in 
d\ — in spite of particle one apparently being a spectator in this interaction. Again, the 
simplest choice 

d\(px,P2\ Q) = id c l (p 1 )d 2 (p 2 + Q)q 2 l2^d 2 (p 2 ) (4.41) 
does not obey the Ward-Takahashi identity. However, if we construct d c according to 

% r 

d^(pi,P2] Q) = - [(dl(pi) + rfi(pl))rf 2 (P2 + Q)q2i2^d 2 (p 2 ) 

+ dl(p 1 )j 2 ,44d c 1 (pi)(d 2 (p2 + Q) + d 2 (p 2 ))} (4.42) 
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then it does obey the WTI (|4.40|) . This can be shown as follows. Contracting the right-hand 
side of Eq. Q4.42|) with Q M gives 



i r 
2 



But, 



(d?(pi) + d\(p 1 ))q 2 (d 2 {p 2 ) - d 2 (p 2 + Q)) 

+dl{ Pl )q 2 2{4' - <) 7 o^( Pl )(d 2 (p 2 + Q) + d 2 (p 2 )) 

2(«S / -KS)7io = ^(pi)- 1 -d?(Pi)- 1 , 



(4.43) 



and so, the Green's function defined by Eq. Q4.42 ) does indeed satisfy the WTI ( |4.40 ). 

It is now possible to combine the results for jj, 2 ^ and d°„ in order to construct the 4D 
current Qq ^ = Gq^ + G Cim corresponding to the free Green's function ( |2.24j ). In this case the 
relevant WTI (again, only writing explicitly the parts where the photon couples to particle 
two) is: 



Q^QIM^ Q) = 2<^i(Pi) + dt( P i))(d 2 (p2) + d c 2 {p 2 )) 

-~<fe(di(pi) + dl{ Pl )){d 2 {p 2 + Q) + dl(p 2 + Q)) + (1 <- 2). 



(4.44) 



Using the forms already constructed for and d c ^ we discover that the 4D current in the 
formalism involving Gc is: 

GoAPi,P2;Q) = 



- [(2di(pi) + ^(pi) + dl( Pl ))(d 2 (p 2 + Q)q 2 ^d 2 {p 2 ) + d\{p 2 + Q)j^d c 2 (p 2 )) 

+ dl( Pl )j 2 ^hidi(Pi)(d 2 (p 2 ) + d 2 (p 2 + Q) + d c 2 (p 2 ) + d\{p 2 + Q))\ + (1 <- 2), (4.45) 

and it obeys Eq. (fj~Hp . 



D. Reduction to 3D and the ET current 

Having constructed a 4D current for the formalism involving Gc that obeys the required 
Ward-Takahashi identity, we can apply the reduction formalism of Section IV B to obtain 
the currents corresponding to the 3D reduction of this 4D theory. To this end, we first note 
that the 4D vertex functions are related to 3D ones, Tet, by Eqs. (|1.15|) and (|1.17|) . Using 
these relations to rewrite the 4D current ( |1.21| ) in terms of the 3D vertex functions shows 
that 



A ETt , = f 1 , ET (P')(^ L ( J P / )^^( P )) r i,ET( J P) 



+ f liET (p')(^(P / )^^(^))ri 1 E T (P). 



(4.46) 



Here, denotes the interaction current obtained by coupling the photon in all possible ways 
within the reduced interaction U . In the same manner as discussed for the 3D version of 
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the Bethe-Salpeter current, it is possible to maintain gauge invariance if the 3D interaction 
U\ is truncated at some definite order in the coupling constant, while the expansions for 
VL l Qq ij VL r and QU^Q similarly are truncated at the same order in the coupling constant. 

Because we truncate our effective interaction at second order, the corresponding con- 
served current in 3D form is obtained by use offi R(2) = l + [U {2) -Ui 2) }g and a corresponding 
expansion for in Eq. (^4.46|) . The result for the 3D current, 

4 2) = f w(P , )(^)ri,E T (P) 

+ r 1 , ET (p')(^(P / )(^ {2) (^ / ) - u^iP'Wl^MP) 

+ f iMp')(G1^ 2 \p) - f/ 1 (2) (P))^o(P))r 1 ,ET(P) 

+ r liET (p')(^o( J P , )^ 7{2) ^o( J P))r 1 ,E T ( J P), (4.47) 

obeys the Ward-Takahashi identity ( f4.15| ). 

Note that = and = are correct at second order in expansions in the 

coupling constant. Consequently if we are considering one-boson exchange interactions the 
only contributions to that are necessary for gauge invariance come from the photon 

coupling to isovector exchange particles and from contact terms due to derivative couplings 
of the mesons to the nucleon. Since these give rise to isovector structures their contribution 
to electromagnetic scattering off the deuteron is, in fact, zero. 



E. Impulse-approximation current based on the instant approximation to ET 

formalism 

Just as in the case of the Bethe-Salpeter equation, if the instant approximation is used 
to obtain from Eq. (|2.25|) the three-dimensional equation that defines the vertex function, 
then a corresponding gauge-invariant impulse current can be constructed by replacing U 
by Ui nst everywhere in the expression (|4.47|) . Assuming that the interaction in the four- 
dimensional equation does not depend on the total momentum, this leads to a particularly 
simple conserved current: 

>W = Finst(£ 7 ,M> r inst- (4.48) 



At this point it appears that we must integrate the rather complicated formula ( f4.45| ) 
over the zeroth component of relative momentum in order to calculate (Qq^)- However, 



the result ( 4.45 ) for Qq„ was constructed in order to obey Ward-Takahashi identities in 



the full four-dimensional theory. It is not, in fact, necessary to use the full result if we 
are only concerned with maintaining WTIs at the three-dimensional level in the instant 
approximation. As remarked following Eq. ( |1 . 1 2| ) , the simpler propagator (di(d 2 + d^)) 
provides the same result as (Go)- Thus, we may construct the corresponding current 

£i 7 nst,>i, p 2 ; P, Q) = i(di{ P i)d 2 (p 2 + Q)j^d 2 ( P2 ) + d l ( Pl )dl{ P2 + Q)j®d°{p 2 )} + (l <- 2). 

(4.49) 

When contracted with Q M this gives 
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Q^im,n(PuP2] Q) = iq2[(di(pi)d 2 (p 2 )) - (d 1 (p 1 )d 2 (p2 + Q)) 

+ {d l {p l )d c 2 {p 2 )) - (d l {p 1 )dl{p 2 + Q))} + (1 <-> 2) 

= <&[<&> (pi, p 2 ; P°) - (SoXpi, p 2 + Q; P° + g )] 
+(1 <->2). 



(4.50) 



(4.51) 



Consequently if a vertex function r inst is constructed to be a solution to Eq. ( |2.26|) then the 
current matrix element defined by 



A 



inst,/^ 



d 3 p 

J2^y 



r inst (p - Q/2; p')g£L(pi p ; Q)r inst (p; p) + (i <-> 2) 



(4.52) 



obeys the current conservation condition ( 4.18Q . The current £/^ stA( is simpler than the full 
ET current and omits only effects stemming from retardation in the current. Since our 
present calculations are designed to provide an assessment of the role of negative-energy 
states and retardation effects in the vertex functions, and these effects stemming from re- 
tardation in the current are expected to be minor, G? nst ^ is used in these calculations. 

V. ELECTRON-DEUTERON SCATTERING: IMPULSE APPROXIMATION 

CALCULATIONS 

Now we are in a position to calculate the deuteron electromagnetic form factors A and B, 
and the tensor polarization T 20 . These are related to the charge, quadrupole, and magnetic 
form factors of the deuteron, F c , Fq, and F M , by the following formulae, 



where 



x 



V 



A-- 
B -- 

^20 = 
^F Q 

3F C ' 
2r] 
~3~ 



p2 , °„2i72 , z '„ r 2 



-V2 



x(x + 2) + y/2 
l + 2(x 2 + y) ' 



Q 



+ (1 + r/)tan 2 



F c 



4M£' 



(5.1) 
(5.2) 
(5.3) 

(5.4) 
(5.5) 
(5.6) 



where Q 2 is the square of the four-momentum transfer. For elastic electron-deuteron scat- 
tering, Q is space-like and rj is positive. In this work all plots and data are quoted at an 
electron angle of 9 e = 70 degrees. 

Deuteron form factors Fc, Fq, and Fm are related to the matrix elements of the current 
A^ discussed in the previous section, taken between the three different magnetic quantum 
number states \M) = | + 1), |0), and | — 1) of the deuteron as follows: 
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F c 
Fq 
Fm 



1 



1 



2r)y/l + rje 
-1 

2?7(1 + r])e 



((0|^°|0) + 2(+l|^°| + l)), 
((0|^°|0}-(+l|^ |+l)), 



i+l\A + \0). 



(5.7) 
(5.8) 
(5.9) 



Matrix elements of A^ are calculated in the Breit frame, with kinematics as shown in 
Fig. [T2|. In order to calculate matrix elements such as those defined by Eq. (|4.52| ) we require 
the vertex functions, or equivalently the wave functions, in the frame where the total four- 
momentum is P = 



' + Q 2 /4, Q/2). However, these are precisely the wave functions 
calculated in Section |T|. Thus, we now take the wave functions constructed for the five 
different interactions of Section [III] at P 2 ranging from to 25 fm -2 and insert them into 
the expression ( |4.52j ). The explicit form of the three-dimensional current Q\ nst ^ is presented 
in Appendix ||. In using any of the interactions obtained with only positive-energy state 
propagation we of course drop all pieces of the operator Q^ nst ^ in negative-energy state 
sectors. 

The general form of a matrix element of the operator A^, using the deuteron vertex- 
function decomposition developed in Ref. |jT9|, is, for coupling only to particle one, 



{M'\Ap\M) 



d 3 p rT ^p' lS ' lP2S2 ,M' 



(27T) 



p + Q/2, P + Q)]*<£ > p ' lP2 (P + Q/2; P + Q) 



J* 



S 1 p±Sip2S2 



(Sb)^(p;P)rf 



P1S1P2S 2 ,M 



with: 



(g y^( p - P )= PlP2T 



and 



p' 1 s[plSlp2S2 



(pi + p2)P° -e 1 — e 2 + \{pi - p 2 )(ei - e 2 ) 



P2u p, ^p\v\)j(:r uPiS1 ^ 



p;P) + (!<-> 2), (5.10) 



(5.11) 



(5.12) 



where the form of Jf i 1 A f lP2 is rho-spin dependent and is given in Appendix |B|. 

The single-nucleon currents used in these calculations are the usual one for extended 
nucleons, defined for n = 1 or 2 by 



(n),p 



qn ( Y n F[ n) m + ^r< U Q v Fr{Q 



(n) t 



2M 



N 



(5.13) 



The appearance of the factor F\ in the first term here implies that the WTI Q^jp —Q is n °t 
satisfied. On the other hand, as pointed out by Gross and Riska a current which obeys the 
original WTI is 



•LL 

J\VT 



7 P + 



2M 



-cr^Q v F 2 (Q 2 )+ F 1 (Q 2 



1 



N 



r 



Q> 



Q 2 



(5.14) 
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The difference between j M and is proportional to Q^. Therefore, for elastic electron- 
deuteron scattering in the Breit frame only j 3 will be affected if we adopt the WTI-obeying 
current 2wt instead of the original extended- nucleon current Recall that in the Breit 
frame j 3 is formally constrained by current conservation to be proportional to j°, and is never 
evaluated. Hence, it transpires that by formally adding a piece to the one-body current j^ 1 
we can obtain a one-body current for an extended object which still satisfies Q^j^ =jQ, 
and yet, in this calculation, leads to the same numerical results as Eq. ( |5.13j ). We choose to 



parametrize the single-nucleon form factors Fi and F 2 via the 1976 Holder fits jEy. Choosing 
different single-nucleon form factors does not affect our qualitative conclusions, although it 
has some impact on our quantitative results for A, B, and T 20 . 

Using such a one-body current the calculation described above conserves the electromag- 
netic current if the vertex function T\ is obtained using an instant potential. However, in 
all other circumstances it violates the Ward-Takahashi identities by omission of pieces that 
are required because of the inclusion of retardation effects in the calculation. This violation 
of the WTI's comes from two sources. 

Firstly, in, for instance, the TOPT calculation, the pieces of the current coming from 
terms in Eq. (|4.13| ) of the form 

(G (P')(^ (2) - k{ 2) )G (P)J,G (P)) 

are not included if Eq. ( |4.52| ) is applied. The inclusion of these "in-flight" contributions 
to the current in our calculations, and the consequent restoration of current conservation 
in the TOPT calculation, will be the subject of a future paper. There these contributions, 
and meson-exchange current effects such as the pn^y and uaj MECs will be calculated. Of 
course, similar "in-flight" pieces are missing from the "retarded ET" calculation. 

Secondly, even if these "in-flight" pieces are included in the retarded ET calculation the 
WTI's will not hold. Pieces will still be missing from the current ( 4.52)) . Specifically, we need 



to use the more complicated form of Qq „ given in Eq. ( |4.45| ). However, calculations where 



(Gq ) is omitted from the instant current matrix element suggest that the total contribution 
of (Gq ) to the observables is small. Since (G^„) itself makes a small contribution to 
observables, using a form of it that only approximately satisfies current conservation is 
expected to have an even smaller effect on our numerical results. 



VI. CONCLUSION 

The ET formalism developed here provides a systematic three-dimensional theory of elec- 
tromagnetic interactions involving relativistic bound states. This is achieved by integrating 
over time components of relative momenta. For the propagator of the theory, this produces 
a form corresponding to zero relative time of the two particles. In order for this formalism 
to incorporate the Z-graphs that are expected in a quantum field theory, it is necessary for 
the propagator to include terms that come from crossed Feynman graphs. The predominant 
terms that are needed have been derived using a form of the eikonal approximation. This 
leads to the ET propagator. In this paper, we discuss a refined version of the ET propagator 



involving a different choice of k 2 from that of Refs. |HI[II|. This new choice of k 2 avoids 



the unphysical singularities which otherwise occur in the three-dimensional ET interaction 
when it is evaluated in a frame where the bound state has large three-momentum. 
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Given a suitable choice for the ET propagator, one may calculate the interaction and the 
electromagnetic currents that must be used with it. A full accounting of both the couplings 
to negative-energy states and the role of retardations in the interaction is thereby obtained. 

When the electromagnetic currents are constructed in our 3D theory the necessary Ward- 
Takahashi identities are clearly satisfied provided that one calculates both the the interaction 
and the electromagnetic current to all orders. However, this is neither practical, nor, if an 
effective hadronic Lagrangian is being used, desirable. A truncation of both the interaction 
and the electromagnetic current is needed. It is not a priori evident that one may truncate 
the expansion for the both the interaction and the electromagnetic current in a way that 
systematically maintains current conservation. We show that this can in fact be done in the 
ET formalism. 

Violations of Lorentz invariance occur when the interaction is truncated. For the range 
of momenta that arise in electron deuteron scattering at TJNAF Q 2 these violations may 
be compensated for by the use of a renormalized interaction. It is found that the required 
renormalization typically is no more than a few percent. Based upon this, we do not expect 
this violation of Lorentz invariance to have a significant effect on results for electron-deuteron 
scattering. 

Calculations have been performed for the impulse approximation, in which we use an 
instant approximation for the electromagnetic current. This current satisfies current conser- 
vation when used with deuteron vertex functions that are obtained with instant interactions. 
We also have used this simple current with vertex functions that include the full retardations 
obtained in the ET formalism. 

Impulse approximation results fall systematically below experimental data for the form 
factors A and B for Q > 3-4 fm -1 . This deficiency of the theoretical calculations at higher 
Q indicates that additional mechanisms beyond the impulse approximation should be sig- 
nificant. However, the existing tensor polarization data are reasonably well described, and 
this is consistent with previous analyses that have shown T20 to be less sensitive to two-body 
currents. 

The role of negative-energy states is found to be not very large. Our impulse- 
approximation numerical results are in closer agreement with those of Hummel and 
Tjon |~5|-[l7] than with those based upon the spectator formalism [Q. Because the ET 
formalism incorporates the relevant Z-graphs in a preferable way, we are confident that 
Z-graphs play a minor role in calculations that are based upon standard boson-exchange 
models of the NN interaction. The role of retardation corrections in the deuteron vertex 
functions also is rather minor. 

Further calculations are needed in order to incorporate the full ET current and the 
meson-exchange currents. 
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APPENDIX A: LEADING EFFECTS IN 1/M OF TIME-ORDERED Z-GRAPHS 



Here we examine the fourth-order graphs for the NN t-matrix in the equal-time for- 
malism and show that the contributions of order 1/M are correctly reproduced in the ET 
equation with instant interactions. 

The fourth-order piece of the equal-time Green's function is: 

(G K^G K^G ) + (G K^G ) (Al) 

where K is the full Bethe-Salpeter kernel. If we consider only positive-energy states on the 
external lines, and decompose the internal nucleon lines according to the different rho-spins 
which are possible, we produce forty-eight graphs. However, if we restrict our attention to 
graphs which contribute at order 1/M or above only the twelve graphs shown in Fig. [13| are 
relevant. Of these, those on the first and third lines arise from the first, iterative, part of 
expression (|A1[) while those on the second and fourth line come from the crossed-box piece 
of Eq. flATp. 

The external Green's functions are now amputated, and the static limit of 

tH^ ++ = (G Q )lX [(GoK^GoK^Go) + (G K^G )] ++ ^ ++ (G >^ (A2) 

is taken. In this limit the first two graphs in Fig. O are infinitely enhanced, as they come 
from the iteration of the lowest-order interaction . Meanwhile the next four graphs give 
energy denominators: 

111111111 . , 

2 ■ , (A3) 



U)\ U)\ + UJ 2 Co>2 ^2 ^1 + ^2 ^1 W\ &l + ^2 W\ 

where U\ {uj 2 ) is the energy of the first (second) meson emitted from nucleon one. These 
graphs do not vanish in the limit M — * oo. They correspond to a shorter-range interaction 
than k[ 2 ^ , and are formally part of the fourth-order three-dimensional kernel . Their ef- 

(2) 

fects are usually subsumed into K\ through the fitting of the parameters in that interaction 
to the NN scattering data. 

The contribution of the remaining six graphs begins at order 1/M. At that order they 
give energy denominators 

111 1 1111 1 



uj\ 2M U2 2M u)\ + UO2 UJ2 k?i u>i + u>2 2M 

111 11 11 11 
LO2 2M iO\ 2M lo\ + U02 uo\ LO2 lo\ + LO2 2M 

Including the appropriate spinors and factors to obtain the full expression shows that the 
sum of these graphs yields: 

3 %+y 1 (i) tt r4^i(2)4%^(iW4^(2)4 \-- 2 ) (~) , (as) 



where Vi(j) is the vertex for the interaction of meson i with nucleon j. Note that in writing 
this expression we have taken some liberties with the ordering of these vertices. However, 
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the error that this produces is of higher order in 1/M. To order 1/M Eq. ( |A5[) reproduces 

the +H > ++ matrix element of the iteration of the instant interaction using the — h sector 

of the ET propagator: 

v + v + K {2) Al A ^ K {2) v + v + (AfTl 

U l U 2 -^inst _2M inst 1 2 ' ^ ' 

Similarly, the 1 •*->• 2 version of expression (|A6| ) gives the six fourth-order particle- two Z- 
graphs correctly to order 1/M. 

Thus we see that using the ET propagator with an instant interaction correctly repro- 
duces the 1/M pieces of the fourth-order Z-graphs in the scattering series. By contrast, the 
Gross or BSLT propagator with an instant interaction would yield a fourth-order Z-graph 
contribution which is only half the correct size. 



APPENDIX B: EXPLICIT FORM OF THE IMPULSE APPROXIMATION 

CURRENT 



In this appendix we present the explicit form of the impulse approximation current which 
yields a current- conserving calculation if vertices obtained using an instant interaction are 
considered. Recall that 



^i 7 nst, At (Pl ) P2 



P, Q) = iidifa + Q)j^d 1 ( Pl )d 2 ( P2 ) + dl( Pl + Q)jVdt( P i)d 2 (p 2 )) + (1 <- 2). 

(Bl) 



Inserting the forms of the propagators and performing the integrals over relative momenta 
we arrive at a form which is most easily expressed in terms of the quantity Jjl lPlP2 , as written 
in Eq. (|5.10| ). In rho-spin conserving sectors it is: 



J p 



3 p. iX &piP2)jpi 



(B2) 



where j M is the object defined for nucleon two in Eq. ( |4.38| ). In our calculation the defined 
there is multiplied by Fi(Q 2 ) in order to take into account the extended-nucleon structure. 
Likewise, j M is the usual extended-nucleon current ( [5.131 ). 
Meanwhile, in rho-spin changing sectors we have: 



J p 



P°- ei -6 2 



Q° + ei + ei ' Q° + ei + ei ' Q° + 2(«?-<) 
P°-e 1 -e 2 

i 3p 



ei - ei 



g + 2(K?-O-e 1 -e' 1 - 



(B3) 



T+- + 



1 + 2 



e 2 



e'i - e 2 



P 



o' 



ei - e 2 



+ 



Q°- 

pOI 



ex -ei Q° + 2(«? - <) + e t + e x > Q°-e 1 - e[ _ 
- e'i - e 2 



3 (j, 



Qo + 2(«§ - <) + ei + ei 



(B4) 



J" 



1 + 2 



P 0/ + ei + e 2 



P 0/ + ej + e 2 

Q° + e x + ei Q0 + 2(«? - K?') - e x - ei Q° + ei + ei 



J// 
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+ 



P ' + e[ + e 2 



go + 2 ( K o _ „o/) _ £i _ e / 



(B5) 



J;: 



P° + ei + e 2 



P + e! + e 2 



P° + ex + e 2 



g°-ei-ei Q° + 2(«?-«S')+ei + ei 



J A 4 



QO + 2(«? - <) + ex + e; 



(B6) 



These reduce to the formulae of Devine |T^] if the appropriate limits are taken. Note that 
in ++ states the current J M is just the usual single-nucleon current (|5 . 1 3 ) . 
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TABLES 

TABLE I. Meson quantum numbers, masses, cutoffs, and couplings as taken from the Bonn-B 
model. Note that the number in brackets in the p row is the p-meson tensor coupling. 



Meson 




T 


Mass (MeV) 


Cutoff (MeV) 


g 2 /^ 


7T 


o- 


1 


138.03 


1200 


14.6 


V 


1- 





548.8 


1500 


5.0 


p 


1+ 


1 


769.0 


1300 


0.95 (6.1) 




1+ 





782.6 


1500 


20.0 


5 


0+ 


1 


983.0 


1500 


3.1155 


a 


0+ 





550.0 


2000 





TABLE II. Sigma coupling required to produce the correct deuteron binding energy in the five 
different models under consideration here. 



Interaction 


States included 


«£/4tt 


Instant 


++ 


8.08 


Klein 


++ 


9.64 


Retarded ET 


++ 


8.39 


Instant ET 


All 


8.55 


Retarded ET 


All 


8.44 
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FIGURES 




FIG. 1. One example of a Z-graph which is not included in the ladder Bethe-Salpeter equation 
scattering series 




FIG. 2. One example of a two-body current that is required in our formalism in order to 
maintain gauge invariance. 




0.0 25.0 50.0 75.0 100.0 



Q 2 (fm 2 ) 

FIG. 3. The form factor A(Q 2 ) for the deuteron. The dash-dotted line represents a calculation 
using a vertex function generated using the instant interaction. Meanwhile the solid line is the result 
obtained with the retarded ET vertex function. In both cases both positive and negative-energy 
sectors are included. 
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FIG. 6. The form factor A(Q 2 ) for the deuteron, as predicted by three different NN models. 
The dotted curve is the instant calculation, the dash-dotted curve is the calculation using the TOPT 
interaction, and the solid curve is the result using a retarded ET interaction. All calculations are 
done with the contributions of negative-energy states dropped. 



q 

CD 



10" 



10' 



10 



10' 



10" 



10' 





Instant calculation: ++ only 
TOPT calculation: ++ only 
Retarded ET calculation: ++ only 












IE 










r J S 













0.0 



25.0 



50.0 

Q 2 (fm 2 ) 



75.0 



100.0 



FIG. 7. The form factor B(Q 2 ) for the deuteron, legend as in Fig. ||. Again, all calculations 
are done in the +- 1- only approximation. 
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1.0 




-1.5 1 — 1 — ■— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
0.0 25.0 50.0 75.0 100.0 

Q 2 (fm 2 ) 



FIG. 8. The tensor polarization T2o(Q 2 ) in electron-deuteron scattering, legend as in Fig 
All calculations are done using only ++ states. 

k', p - P 2 k, 



FIG. 9. The crossed-box graph, showing the momentum labels used in the text. 



FIG. 10. The two time-ordered perturbation theory graphs for one-pion exchange. 
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1.25 




0.95 ^ — ' — — ' — 1 

0.0 10.0 20.0 30.0 40.0 50.0 

P (fm ) 



FIG. 11. A plot showing the eigenvalue, A, of the homogeneous integral equation (|2.26| ), as a 
function of the total three-momentum of the two-body system, for different choices of the interaction 
U\. The dot-dashed (dotted) line is the result when U\ is chosen to be an instant interaction and 
all (only +- 1-) states are included in the calculation. The solid (long dashed) line is the result 
when U\ is chosen to be the retarded interaction defined by Eq. ( |2.33| ) and all (only +- 1-) states are 
included. Finally, the short-dashed line is the result from the TOPT interaction, calculated with 
+- 1- states only. 
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(0,Q) 



(E.Q/2) 




(E.-Q/2) 



FIG. 12. The kinematics of electron-deuteron scattering, for a momentum transfer Q = (0, Q), 
in the Breit frame. 



+ (1 and 2 interchanged) 

FIG. 13. The graphs which contribute to the three-dimensional NN t-matrix at fourth order in 
the coupling, up to order 1/M. The graphs on the first and third line come from the NN Feynman 
box graph and those on the second and fourth line come from the NN Feynman crossed-box graph. 
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